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ABSTRACT 

Interacting double white dwarf binaries can give rise to a wide variety of astrophysical outcomes 
ranging from faint thermonuclear and Type la supernovae to the formation of neutron stars and stably 
accreting AM Canum Venaticorum systems. One key factor affecting the final outcome is whether 
mass transfer remains dynamically stable or instead diverges, leading to the tidal disruption of the 
donor and the merger of the binary. It is typically thought that for low ratios of the donor mass to the 
accretor mass, mass transfer remains stable, especially if accretion occurs via a disk. In this Letter, 
we examine low mass ratio double white dwarf binaries and find that the initial phase of hydrogen- 
rich mass transfer leads to a classical nova-like outburst on the accretor. Dynamical friction within 
the expanding nova shell shrinks the orbit and causes the mass transfer rate to increase dramatically 
above the accretor’s Eddington limit, possibly resulting in a binary merger. If the binary survives 
the first hydrogen-rich nova outbursts, dynamical friction within the subsequent helium-powered nova 
shells pushes the system even more strongly towards merger. While further calculations are necessary 
to confirm this outcome for the entire range of binaries previously thought to be dynamically stable, 
it appears likely that most, if not all, interacting double white dwarf binaries will merge during the 
course of their evolution. 

Subject headings: binaries: close— novae, cataclysmic variables— nuclear reactions, nucleosynthesis, 
abundances— supernovae: general— white dwarfs 


1. INTRODUCTION 

Double white dwarf (WD) binaries are formed via com¬ 
mon envelope evolution in binary systems. The emission 
of gravitational waves in the closest double WD binaries 
shrinks their orbits and causes mass transfer between the 
WDs within a Hubble timescale. The resulting interac¬ 
tion leads to a broad range of outcomes including ex¬ 
plosive transient phenomena and the formation of exotic 
classes of stars. 

Previous work has found that the dynamical stability 
of mass transfer in these binaries, and thus their ulti¬ 
mate fates, depends on the mass ratio of the donor mass 
to the accretor mass, q = M 2 /M 1 , and on the timescale 
for an gular momentum feedback from the accret or to the 
orbit (Marsh et al. 20041 Kremer et al. 2015). How¬ 
ever, systems are generally assumed to be dynamically 
stable if accretion occurs via a disk, as torques on the 
extended lever arm of the disk should yield efficient an¬ 
gular momentum feedback. In this Letter, we challenge 
this assumption by showing that even disk-accreting dou¬ 
ble WD binaries with extremely low-mass He WD donors 
will merge unstably. 

The overlying non-degenerate H-rich layers surround¬ 
ing these low-mass WDs alters the initial phase of mass 
transfer compared to s implified pure He W D donors 


(D’Antona et al. 
to classical nova-. 


2006 Kaplan et al. 2012) and leads 
ike 011 tbursts (IShen et al.||2018|) . Sub- 
:lyl 


sequent accretion of the underlying He-r ich WD core 
yields analogous He-powered nova episodes (Shen & Bild- 
sten 2009). Here, we consider the impact of these ex¬ 


panding nova shells on the binary orbit. We calculate 
the evolution of two fiducial double WD binaries with 


1 Department of Astronomy and Theoretical Astrophysics 
Center, University of California, Berkeley, CA 94720, USA. 


masses 0.2 + 0.5 and 0.2 + 1.0 M 0 and find that these 
binaries, which were previously assumed to form sta- 
bl y mass transfer ring AM Canum Venaticorum systems 


(Kilic et al. 


20141, instead merge due to dynamical fric- 


tion within the accretor’s expanded outer layers. In Sec- 
tion[2] we estimate the shrinking of the orbital separation 
due to mass ejection via dynamical friction within the 
common envelope formed by the nova shell. In Section 
[3j we provide a prescription for how the donor’s mass 
transfer rate responds to the decrease in orbital separa¬ 
tion. In Section 4j we describe our fiducial calculations, 
and we consider the implications of our findings in Sec¬ 
tion [5j 


2. COMMON ENVELOPE ENERGETICS 

We begin by estimating the change in the binary sepa¬ 
ration, a, if some or all of the energy to eject the expand¬ 
ing nova shell comes from the orbit via dynamical friction 
of the two WDs with the comm on envelope. This i s at 
best a rough approximation (see Ivanova et al. [2013 and 
references therein); however, we use it here for simplic¬ 
ity and defer analysis with alternative common envelope 
physics to future work. 

The specific energy required to eject material to infin¬ 
ity from the Roche surface of the accreting WD is 
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where the total mass of the system is M t = Mi + M 2 , 
and o;li is the distance from the center of the donor to 
the L\ Lagrange point in units of the binary separation. 
The function f(q) varies from 1.4 — 1.9 for q = 0.1 — 
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0.5. While this specific energy changes slightly as mass 
is being ejected from the system, it will remain correct 
to first order in the ratio of the ejecta mass, M e j, to 
the total mass and thus is sufficient for our purposes of 
estimation. 

Assuming that a fraction oce of the energy to eject M e j 
comes from the binding energy of the orbit, we obtain 
the change in binary separation following the common 
envelope event, 


da 
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which is again correct to first order in (M e j/M t ). Note 
that because of the factor of 2/q, the ratio da/a can be 
more than an order of magnitude larger than 


3. PRESCRIPTIONS FOR THE ACCRETION RATE 

Dynamical friction during the common envelope event 
reduces the binary separation, and thus increases the 
Roche overfill factor, A = R? — Rl 2 , where i? 2 is the 
radius of the donor star, and R \, 2 is the donor’s R oche 
radius, which is approximated by Paczynski (1967) as 
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Figure 1. M2 (top panel) and A (bottom panel) vs. M trans f erre d 
for our fiducial 0.2+ 0.5 and 0.2+ 1.0 Mq binaries, as labeled. Tri¬ 
angles mark when the H mass fraction of the transferred material 
drops below 1%. These time-dependent calculations do not take 
into account the binary’s response to the nova outburst, which oc¬ 
curs at the bullets in each panel, so the evolution after the nova and 
subsequent common envelope is merely shown for completeness. 


As a result of the increased overfill factor, the accretion 
rate rises. In order to estimate the new perturbed accre¬ 
tion rate, we require a relation for M 2 (A). Several such 
prescrip tions exist , including the isothermal approxima¬ 


tion (Ritter 1988); however, as pointed out by Marsh 


et al. (2(J(J4), this approximation breaks down when A is 


larger than the pressure scale height at the donor’s pho¬ 
tosphere, which will indeed be the case for our post-nova 
binaries. 

We therefore use the adiabatic M prescription 


(Paczynski & Sienkiewiczl 1972; Savonije 1978 

Webbink 

1984 Marsh et al. 

2004), which assumes an equation ot 


tential surfaces. However, because we are interested in 
the early phase of mass transfer involving the donor’s 
outer, non-degenerate layers, we cannot assume a con¬ 
stant polytropic index for all donor material as the pre¬ 
viously cited studies do. Instead, we follow the more 


general expression of Ge et al. (2010), which allows for 
a varying polytropic index within the donor at the com¬ 
plicating expense of requiring a stellar structure integral. 
Their approximation to M 2 (A) is 


M 2 (A) ~ 2t iRlh(q) 
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where c s is the sound speed, the factor h(q) is 


c s pdx( 4) 


4. FIDUCIAL EXAMPLES 

We now utilize the relations in Sections [2] and [3] to 
calculate the evolution of two fiducial double WD bina¬ 
ries consisting of an extremely low-mass 0.2 M 0 He WD 
donor and 0.5 and 1.0 M 0 C/O WD accretors. The WDs 
are construct ed and evolved with the stellar evolution 
code MES+ (Paxton et al.||2011 2013), beginning with 
initial nietamcities of 0.02 and mam sequence masses of 
0.8 and 7 M 0 for the 0.5 and 1.0 M 0 C/O WD mod¬ 
els, respectively. These WDs are taken from the set of 
default WD m odels included with MESA, which employ 
Reimers (1975) and Blocker (1995) winds to remove mass 
on the red giant and asymptotic giant branches, leaving 
hot degenerate cores. We begin with these cores, turn 
on chemical diffusion, and allow them to cool for 10 10 yr 
to approximate the long main sequence timescale of the 
low-mass companion that will give rise to the 0.2 M 0 He 
WD. 

We construct our fiducial 0.2 M 0 WD by halting 
the evolution of a 1 M 0 star once it evolves off the 
main sequence and expands to a radius of 4 Rq. The 
“cno_extras_ol8_and_ne22” nuclear network is utilized, 
and chemical diffusion is active throughout the simula¬ 
tion. Mass is removed until only 0.2 M 0 remains, and 
the resulting hot WD is allowed to cool for 10 9 yr. This 
yields a He-ricli core mass of 0.194 M 0 and a H-rich shell 
of 0.006 M e . 
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4.1. Deriving the evolution of the mass transfer rate 

In order to calculate the mass transfer history of these 
two binaries, we first ignore the accretor’s response to 
accretion and focus only on the donor’s response to 
mass removal. Since mass transfer is not extremely 
rapid, it is crucial that we account for the donor’s non- 


and we have assumed A <c i? 2 - The function h(q) ranges 2 http://mesa.sourceforge.net, version 7184; default options used 

from 17—14 for q = 01 _05 unless otherwise noted 
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adiabatic thermal response. We first remove mass from 
the 10 9 -year-old 0.2 M 0 donor at a constant rate of 
1CU 9 Mq yr -1 for the donor’s entire history, then cal¬ 
culate a new mass transfer history based on the donor’s 
radial response. This is done assuming orbital angular 
momentum loss is due to gravitational waves and that the 
angular momentum lost from the donor is transferred ef¬ 
ficiently back to the orbit. This limit is likely applicable 
to the 0.2 + 1.0 Mq binary b ecause accretion occurs via 
a disk (Lubow & Shu| 1975) but is only an approxima¬ 
tion in the 0.2 + 0.5 Mq system, which undergoes direct 
impact accretion. These assumptions yield an accretion 
rate of 


M2 

M 2 


32G 3 M\M 2 M t (5 , £ 
5 c 5 a 4 


6 + 2 q 
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where £ = d In R 2 /d\n M 2 . We then restart the donor 
evolution with the new time-dependent M 2 , use the 
donor’s response, including the changing value of £, to 
calculate a second time-dependent M 2 , and iterate once 
more. The maximum difference in M 2 histories between 
the penultimate and final iterations is < 1% after the 
initial ~ 10 -4 Mq has been transferred. The converged 
M 2 histories are shown in the top panel of Figure [l] as a 
function of the transferred mass, M trans f e rred , for our two 
fiducial binaries. 

Also shown in the bottom panel of Figure [T] is the 
Roche overfill factor corresponding to M 2 at a given 
Mt rans ferred) as defined in equation Q. It ranges from 
~ 10 5 cm while the non-degenerate outer layers of the 
donor are transferred to -C 10 5 cm as the degenerate He 
core is excavated. 


4.2. The binary’s response to the expanding accretor 

The M 2 history, including the changing accreted com¬ 
position, is then used to evolve the accreting WDs. 
Chemical diffusion is active in these calculations as be¬ 
fore. After 1.8 x 10 -3 M 0 has been transferred to the 
0.5 M 0 accretor, nuclear burning triggers convection in 
the outer 2.8 x 10 -3 Mq, and a classical nova-like out¬ 
burst p ower ed by H-burning begins. As we will discuss in 
Section|4.3l the 0.2 + 1.0 Mq binary will likely remain dy¬ 
namically stable after its first H-powered outbursts. We 
thus instead begin the accretor’s evolution for this binary 
after the initial H-rich 0.007 Mq has already been trans¬ 
ferred and the accreted material transitions to > 0.99 
He by mass. Convection then begins on the 1.0 Mq ac¬ 
cretor after an additional 3.3 x 10 -3 M 0 of He-rich ma¬ 
terial is accreted, resulting in an expanding envelope of 
2.9 x 10 -3 M 0 . Bullets and squares in Figure[l]show the 
accretion rates and Roche overflow factors for the first 
H-powered and He-powered novae, respectively. 

As the convective layer heats and expands, it overflows 
the accretor’s Roche lobe and forms a common envelope 
surrounding the binary. Even though the majority of 
the convective layer is expelled, the ejection of even a 
small fraction of the shell via dynamical friction is suf¬ 
ficient to increase the donor’s Roche overfill factor and 
consequently M 2 by a tremendous amount. 

Figure [ 2 ] shows the donor’s M 2 (A) relation at the mo¬ 
ment when convection is initiated on the accretor. Dia¬ 
monds demarcate M 2 (Aq) just prior to the common en¬ 



Figure 2. M 2 vs. A for our fiducial binaries, as labeled, at the 
moment when H-burning and He-burning on the 0.5 and 1.0 Mq 
accretors, respectively, initiates convection. Diamonds mark the 
state of the donor prior to the ensuing common envelope event. 
Triangles and squares denote the resulting A and M 2 if orbital 
energy is used to eject 10% and 100% of the ejecta, respectively. 
Crosses represent the Eddington accretion rates for the 0.5 and 
1.0 Mq accretors. 

velope event. The overfill factors resulting from the de¬ 
creased orbital separation, assuming = 0.1 and 1.0, 
and the corresponding mass transfer rates are shown as 
triangles and squares, respectively. It is clear that the 
ejection of a relatively small envelope via a common en¬ 
velope has a drastic effect on these systems, resulting 
in accretion rates that are highly super-Eddington with 
respect to the accretors; the Eddington accretion rates, 
-Mecm) are shown as crosses. 


4.3. The post-common envelope evolution 

We now calculate the evolution of our fiducial systems 
in their post-common envelope, super-Eddington accre¬ 
tion states. To simplify matters, we approximate the 
donor’s M 2 (A) relation as a power-law: M 2 oc A 6 . A 
value of b = 3 is app ropriate fo r pure ly degenerate ob¬ 
jects, as utilized by|Marsh et ah (200 1|) and others. Even 
though the donors in our fiducial binaries are not fully 
degenerate, the results in Figure [2] imply that b ~ 3 in 
these cases as well. With this simplification, the evolu¬ 
tion of M 2 can be written as 


dM 2 / dt bR 2 
Mo ~ A 


rapid 


1 \ M 2 d 1 M t 
3 ) W 2 ~ a + 3 M~ t 


■ ( 8 ) 


We have assumed that A <C R 2 and used the |Paczynski| 
(1967) Roche lobe relation. The donor’s radial response 
to the extremely rapid mass loss that follows the common 
envelope event is £ rap id- 

Because the post-common envelope accretion rate is so 
highly super-Eddington, we further assume that M t = 
M 2 . Following a similar energetic argument to that used 
in Section [2j we find 


~ = TT [ 2 q ce (1 + q)f(q) + 1 ], 

a 1VI2 


(9) 
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where f(q) is defined in equation ([Tj) . The evolution of 
M 2 now becomes 


1967); however, since we now expect such systems to 


( iM 2 /dt 
Mi 


bR 2 M 2 
A M 2 


Crapid - 2o:ce(1 + q)f(q) - / 

1 + q 


The radial response of the donor in the 0.2 + 0.5 Mq 
binary during the first nova event is initially £ ra pid ~ 2 
and becomes negative once the electron-degenerate core 
is exposed. The term in brackets is thus negative, and 
the accretion rate becomes larger. However, for the ini¬ 
tial H-rich outbursts of the 0.2+1.0 Mq binary, the donor 
still possesses a relatively large non-degenerate, radiative 
outer layer, and £ rap id 1. Asa result, the donor is able 
to rapidly shrink back within its Roche lobe, and the ac¬ 
cretion rate decreases back to its equilibrium value. This 
system will thus avoid a catastrophic merger during the 
H-rich nova phase. However, after further mass transfer 
exposes the He core, the donor’s radial response to the 
accretor’s subsequent expanding He nova shell becomes 
negative, and the post-nova accretion rate diverges as in 
the 0.2 + 0.5 Mq case. This is the reason we only follow 
the H-poor evolution for the 0.2 + 1.0 Mq binary in the 
previous sections: the system may survive the initial H- 
rich outbursts but will still merge after the first He-rich 
nova. 

With our assumption of a power-law relation for 


M{ A), equation (10) can be integrated to find that these 
systems merge on a timescale that is even shorter than 
the naive post-common envelope accretion timescale, 
M 2 /M 2 , by multiplicative factors of (A/i? 2 ), (1/6), and 
the inverse of the term in brackets, which ranges from —2 
to —6 depending on the value of £ rap id- This is roughly 
the same expression as the timescale for a mildly per¬ 
turbed double WD binar y to retu rn to its e quilibrium 


mass transfer rate (Marsh & Nelemans 2005 Stroeer & 


Nelemans|2009 ), although the accretion rate in our fidu¬ 


cial binaries starts at a much higher value and diverges 
due to the destabilizing effect of dynamical friction. As 
a result, the 0.2 + 0.5 and 0.2 + 1.0 Mq binaries merge 
within hours of the common envelope event. 

5. CONCLUSIONS AND IMPLICATIONS 

In this Letter, we have shown that even disk-accreting 
double WD binaries with extreme mass ratios will merge 
once common envelope effects during the initial phases 
of mass transfer are taken into account. The ejection 
of the accretor’s expanding nova shell causes the binary 
separation to shrink and the donor to dramatically over¬ 
fill its Roche lobe, resulting in highly super-Eddington 
mass transfer rates that lead to a merger. Further pa¬ 
rameter studies covering a range of donor and accretor 
masses and ages will be required to make definitive con¬ 
clusions regarding all double WD binaries; however, it 
appears likely that most, if not all, interacting double 
WD binaries will merge. 

One important consequence of this work concerns the 
contribution of double WD binaries to the population of 
AM Canum Venaticorum (AM CVn) systems, which con¬ 
sist of a possi bly degenerate He donor and a WD accretor 
(see Solheim 2010 for a recent review). One formation 


merge prior to a prolonged phase of He-rich mass trans¬ 
fer, it appears that the AM CVn population must come 
from a different channel, such as the He -burning star 
(Savonije et al.|1986l Iben fe Tutukov 1987) or the H star 


( 10 ) 


channel (Tutukov et al 


■ | 19851 

e doub 


Podsiadlowski et al. 2003). 


I 


The elimination of the double WD formation channel 
would help to resolve the large discrepancy betwe en t he 


predicted and observed AM CVn sp ace densities (Nele¬ 


mans et al. 2001 Carter et al. 2013), which lends some 


support to our findings. 

A related outcome involves the proposed class of “.la” 
supernovae arising from He detonations in stably ac¬ 
creted He shells during the early lives of double W D bi- 


sten et al. 2007 Shen & Bildsten 

2009, 

Shen et al. 2010) 

Such systems now merge before the d 

onor stably trans 


fers a sufficient amount of He to trigger a detonation in a 
convective He-burning shell. However, He detonations in 
double WD binaries may still occur during the violence 
of double WD mergers if the do nor is a He WD or has 



sequently 

onate in a double detonation Type la supernova (SN la; 
Livne 1990 Shen & Bildsten 2014), the naked He deto¬ 
nation may stiff appear as a .la” supernova. 

The post-nova common envelope events also affect the 
evolution of higher mass ratio double WD binaries that 
would have merged even in the absence of dynamical fric¬ 
tion within the expanding nova shells. In |Shen et al.| 
(20131, we considered a 0.4 + 1.0 Mq binary and exam¬ 
ined the effect of the pre-merger nova episodes on the sur¬ 
rounding interstellar medium. The dense swept-up shells 
act as absorbing screens for a possible subsequent SN la, 
matching recent observation s of blueshifted and/or vari¬ 
able Na abs orption features (Patat et al.|2007 Sternberg 
et al. 2011). However, our calculations did not account 


tor the effect of the common envelopes on the binary evo¬ 
lution. The ejected envelopes in these higher-mass bina¬ 
ries are smaller than for this Letter’s fiducial binaries, so 
the systems may survive the dynamical friction phases, 
but the number of ejected shells before the m erger may 
be lower than calculated in Shen et al. (J2013|). 

We conclude by showing a schematic diagram of in¬ 
teracting double WD binary outcomes in Figure [3j with 
evolutionary phases in chronological order in each do¬ 
main. The long-term evolution of merger remnants that 
do not explode as SNe la during their mergers is as¬ 
sumed to fo llow the basic framework outlined in IShenl 
et al. (2012). He-rich subdwarf B and O stars are ab¬ 
breviated as “sdB/sdO”, R Coronae Borealis stars as “R 
CrB”, iron core-collapses as “Fe CC”, and neutron stars 
as “NS”. Degenerate O/Ne cores that reach high densi¬ 
ties are assumed to undergo runaway electron-captures 
and form a neutron star via accretion-induced collapse 
(AIC), although the standard picture of accretion from a 
companion star does not apply here. SNe la are assumed 
to occur via a core C detonation triggered by a s urface He 
detonation during mergers with Mi > 0.8 M ^ ([ Pakmor 


et al.||2013 Shen fe Bildsten||2014 Shen fe Moore|[2U14| ) 
or via a directly triggered (J deton ation during “violent 
merg ers” of high-mass C/O WDs (Pakmor et al. 2010 


channel involves a low-mass He WD donor (Paczynski 20121 that do not possess sufficiently massive He surface 
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Figure 3. Schematic of interacting double WD binary outcomes. See Section [5] for details. 


layers. 

Boundaries between regions with differing outcomes 
are largely qualitative. Phases in parentheses may not 
occur for all binaries in the region. Dashed lines de¬ 
lineate compositional boundaries between He and C/O 
(~ 0.5 Mq) and C/O and O/Ne (~ 1.1 M 0 ) WDs. 
In the binaries below the dotted line, accretion would 
remain sub-Eddington and they could avoid mergi ng if 


no post-nova common envelopes occurred (Marsh et al. 


2004|. However, as we have shown in this Letter, It Is 


likely that even these binaries merge due to dynamical 
friction within common envelopes. 
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